Future e^e Colliders Sensitivity to Hhh Coupling and CP Violation 
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We perform a complete simulation of the process e'^e" bbvu, where v can be an electron, 
muon or tau neutrino, in the context of a general Higgs coupling to b quarks. We parametrize the 
Hbb coupling as -^(a + i^sb). Taking into account interference effects between pure Higgs and 
Standard Model contributions, we find that sensitivities of the order of 2% and 20% can be obtained 
at a future e"^e~ collider for deviations of the a and 6 parameters respectively from their Standard 
Model values. Combining our analysis with an independent measurement of F^f^jg can provide 
evidence about the CP nature of the Higgs sector. 

PACS Categories: 98.80.Cq 

I. INTRODUCTION 

The origin of fermion masses and mixings is one of 
most important issues in particle physics. In the Stan- 
dard Model (SM), the Higgs field alone is responsible for 
the electroweak symmetry breaking and mass generation. 
The SM, however, is incomplete and a thorough study of 
the coupling of the remnant Higgs boson (or in fact the 
lightest (pseudo)scalar boson) to fermions can provide 
hints on the actual mechanism of mass generation. 

In a recent paper, we have investigated the possibility 
of detecting deviations from the SM in the Higgs cou- 
plings to T-leptons at future e+e" colliders. In this 
letter, we expand our analysis to the case of the Higgs 
couplings to 6-quarks, which has a better potential in 
principle due to the larger Yukawa coupling. 

For definiteness, we concentrate on the determination 
of the (pseudo)scalar-6-6 coupling at a Linear Collider 
with a center-of-mass energy of y/s = 500 GeV and ac- 
cumulated luminosity of 1 ab~^, based on the TESLA 
design Q|. We will assume that this particle has already 
been discovered at the Large Hadron Collider but a de- 
tailed study of its couplings is missing. 

We take into account all relevant contributions to the 
process e~^e~ — > bbvi?, where v can be an electron, muon 
or tau neutrino. In particular, weak gauge boson fusion is 
the dominant contribution to the subset of diagrams con- 
taining the Higgs boson for AIh < 180 GeV at ^/s > 500 
GeV. 

In extensions of the SM with extra scalars and pseu- 
doscalars, the lightest spin-0 particle can be an admixture 
of states without a definite parity. Hence, we parametrize 
the general Hbb coupling as: 
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where v = 246 GeV, nib is the 6— quark mass and a = 1, 



FIG. 1. The total cross section e+e — > bbi^v dependence 
on Aa (solid line) and b (dashed line). 



We will present results considering a and b as indepen- 
dent parameters and also for the cases of fixed a = 1, 
free b and fixed 6 = 0, free a. We will see that there is a 
region of insensivity around circles in the a — b plane since 
we can't at this level of analysis disentangle the effects of 
a and b. 

The cross section for the process e+e" bbi/D is sen- 
sitive to terms proportional to a, which comes from the 
interference with non-Higgs contributions, and and b^ 
from pure Higgs contributions. Therefore, we can search 
for deviations from the SM prediction which could arise, 
for instance, in supersymmetric models. 

The total SM cross section for the process e+e^ 
bbvv is of the order of 180 fb for Mh = 120 GeV and 
at a/s = 500 GeV, being dominated by the Vf. final 
state, because of the many additional diagrams allowed 
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in this case. In particular, only in this channel the weak 
gauge boson fusion diagram is allowed and it provides 
an important contribution. For comparison, the process 
e+e" HVf.Ve is of the order of 100 fb at = 500 
GeV. 

In Fig. |l| we show the dependence of the total cross 
section, summed over the three neutrino species, on the 
parameters Aa = a — 1 and h. The dependence in h is 
symmetric since only terms proportional to contribute 
to the cross section. On the other hand, the dependence 
on Aa is asymmetric due to interference with Standard 
Model processes, which leads to the presence of linear 
terms in Aa. We should point out that the Aa depen- 
dence is more prominent in the bb case than in the pre- 
viously studied r+r^ jl]. Our goal is to see how sen- 
sitive the experiments performed at the next generation 
of e^e" colliders will be in the determination of these 
parameters. 



—aie^e — i hhvv) 



COS Beh COS Bsh 

FIG. 2. The differential cross section over cos 0^6 (in pb) 
for the process e^e~ — > hbv^i7^ and e~ — > }A)vv, summed 
over the three neutrinos. Solid circles are the full SM result 
while crossed circles are the contribution from Higgs only (in- 
cluding interference effects). 



II. ANALYSIS AND RESULTS 



We performed our Monte Carlo simulation by gener- 
ating observables represented as series in the a and b 
couplings multiplied by kinematical factors: 
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where O is any observable and the A^ terms are purely 
kinematical structures which do not contain any a and 
b dependence and results from the amplitude squaring 
and phase space integration. These Ai structures are 
the subject of the Monte Carlo simulation and the a and 
b couplings can be varied without the necessity to re- 
simulating the data for each (a, b) point. In our particu- 
lar case, A^ = Ai = 0. 

The event sample reproducing the expected statistics 
at TESLA was generated using our Monte Carlo pack- 
age while the detector response was simulated with the 
code SIMDET version 3.01 We assume an efficiency for 
6— jet pair reconstruction of Ebb = 56 %, which is based 
on the 5-tag algorithm, as assumed in ref. ||^. In our 
simulations we used Mh — 120 GeV. 

In Figure |^ we show, for comparison purposes, the dif- 
ferential distribution in cos^eb) where Ogb is the scattering 
angle between the 6-jet and initial beam directions, for 
the total SM contribution and for the Higgs contribution 
only (including interference with SM). To illustrate the 
importance of the final state, we plot both the contri- 
bution of only the final state (which is basically the 
same as Vr final state) and the total contribution from 
the three neutrinos. We can see that in the first case the 
Higgs contribution is small but in the total contribution 
it is comparable to the full SM result. 



FIG. 3. Contribution of Higgs diagrams to the differential 
cos^eb distribution in the e"'"e~ — > bhvv' process for Standard 
Model [a = l,h — Q, black dots) and a = 0.5, h = 0.5 (crossed 
dots). 

In order to demonstrate the effect of different values of 
the parameters a and 6, we show in Figure || the cos Oe.b 
distribution arising from the Higgs contribution for the 
SM (a = 1, 6 = 0) compared to the case with a = b = 0.5, 
with only neutrinos and with the three neutrinos. We 
see that the shapes are very similar, as expected, but the 
levels can be noticeable different. 

As for the possible contribution from background pro- 
cesses, like e'^e~ e~^e~ZZ e'^e~bbi>D (with e-pair 
lost), e+e^ -> wW+W- vvbbvv, e+e" ZZZ 
bbvDvv, etc., the cross sections of these processes are ei- 
ther small, or they can be significantly supressed down 
to levels of 0.2 fb g. 

Another important aspect is the assumption about the 
detector performance and possible sources of the system- 
atic uncertaintes. We include the anticipated systematic 
errors of 0.5% in the luminosity measurement, 1% in the 
acceptance determination, 1% in the branching ratios, 
and 1% in the background substraction, and assume the 
Gaussian nature of the sytematics. To place bounds on 
the Hbb couplings, we use a standard x^-criterion to anal- 
yse the events. After various kinematical distributions 
were examined, we found that the most strict bounds 
are achieved from cos 9eb distribution by dividing the dis- 
tribution event samples into 10 bins. The experimental 
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error Acr*^^^ for the i*'' bin is given by: 



III. CONCLUSIONS 



where 



— V °syst + "stat 



(3) 
(4) 



and S'^y^f is the sum in quadrature of the systematic un- 
certainties mentioned above. 

In Fig. ^, we present our final results for a TESLA-like 
environment with a center-of-mass energy of 500 GeV 
and for Mh = 120 GeV. 

We investigated three possible scenarios for the lumi- 
nosities: 100 ih-\ 1 ab~i and 10 ab~^ The allowed 
region for independent Aa and b parameters at 95% con- 
fidence level is the area between the circles. The hori- 
zontal bands are the allowed region for the b parameter 
keeping a = 1 . The vertical bands are the allowed region 
for the Aa parameter keeping 6 = 0. 
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FIG. 4. The allowed regions for Aa and b parameters at 
95% confidence level for £. = 100 fb~^ (long-dashed lines), 1 
ab~^ (solid lines) and 10 ab~^ (short-dashed lines). The area 
between the circles is the allowed region for independent Aa 
and b. The horizontal bands are the allowed region for the 6 
parameter keeping Aa = 0. The vertical bands are the the 
allowed region for the Aa parameter keeping 6 = 0. 
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The bounds that can be obtained at 95% confidence 
level are: 

- 0.041 < Aa < 0.039 for C = 100 fb"!; 
-0.026 < Aa < 0.027 for £ = 1 ab^^; 
-0.024 < Aa < 0.024 for £ = 10 ab^i. 
for the case of 6 = and free Aa and 



0.28 <b< 0.28 for C 
-0.23 <b< 0.23 for C 
-0.22 <b< 0.22 for C 



100 fb-i; 
1 ab-i; 
10 ab-^ 



(6) 



for the case of Aa = and free b. 

These results are up to an order of magnitude bet- 
ter than the limits obtained in a similar manner for 
r— leptons, mainly because of the larger Yukawa coupling 
in case of &— quarks and higher sensitivity of the process. 

These results can be roughly scaled for moderate vari- 
ations in the Higgs boson mass around 120 GeV by mul- 
tiplying the bounds by a factor (M///120 GeV)^. 



We have performed a complete analysis of the sensitiv- 
ity to new Hbb couplings from the process e^e~ — > bbvD 
at the next generation of linear colliders. These new cou- 
plings are predicted by many extensions of the Standard 
Model. We showed that forthcoming experiments will 
be able to probe deviations of Hbb coupling. The weak 
gauge boson fusion process is instrumental for achieving 
such a precision. For a TESLA-like environment, we are 
able to constrain the couplings at the level of a few per- 
cent for the a parameter (for fixed b) and tens of percent 
for the b parameter (for fixed a). These results are com- 
parable to the study performed in , where a global fit 
analysis for C — 500 fb^^ and v^s = 500 GcV has re- 
sulted in a relative accuracy of 2.2% in the gnbb Yukawa 
coupling. 

We would like to comment some aspects of the future 
measurements. Let us assume that the Higgs data an- 
ticipated from the new collider experiments will reveal 
a deviation from the SM predictions. In addition, sup- 
pose that one has an independent measurement of the 
partial width F^^^f, (for instance, from on-mass-shell 
Higgs production in a muonic colhder). It easy to see 
that, in our parametrization, r^^f,^ ~ (a^ -I- &^), while 
the observables we studied have the following parameter 
dependence 
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^ Ao + a- Ai+a^ ■ A2+b^ ■ A3 



Combining our results and those from F^^f,^ one can 
separate a and b contributions and obtain an explicit in- 
dication of CP violation in the Higgs sector. 
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Hi Rogerio, thank you for the info concerning the forthcoming conference. I will apply. 
Now about our paper. Please find below some results for ^/s = 500 GeV and the Higgs 
mass of 120 GeV. 

The cross-sections of the muonic- o tau- neutrino channel is 1.008- 10~^ pb, while total 
is 0.175 pb. 

In figures attached as PS-files one finds: 

Fig. 1. The e+e^ bhvv cross-section dependence on Aa and h parameters. 

Fig. 2a. Differential distribution of the e+e" — > 66/^;i(r)t'/j(T) process over M^j. Solid 
points is the total SM contribution, while crossed points is the contribution from 
Higgs containing diagrams only. 

Fig. 2b. The same as Fig. 2a, but for case of all three types of neutrinos. 

Fig. 3a. Differential distribution of the e+e" bbi'^{T)^^.{T) process over Mjj. Solid 
points is the contribution from Higgs containing diagrams of SM, while crossed 
points is the case of a = 0.5 and b = 0.5. 

Fig. 3b. The same as Fig. 3a, but for case of all three types of neutrinos. 

Fig. 4. Allowed parameters region (area inside the corresponding contour lines) for 95% 
C.L. for the case of: 

1. independent a and b parameters (curled contours); 

2. b = (vertical bands); 

3. a = (horizontal bands). 

Long-dashed, solid, and short-dashed lines represent the cases of 100 fb~^, 1 ab~\ 
and 10 ab~^, correspondingly. 

The efficiency for b-jets reconstruction is 0.75. All other systematics is the same as in 
the tau-tau case. The question is: should we describe the procedure of the beackground 
reduction etc.? We can simply relfer to paper by K. Deusch and N.Meyer, LC-PHSM- 
2001-025. In this paper all the details are perfectly described and the case coincide with 
ours exactly. What do you think? 

Andre 



